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Abstract. We summarize several new developments in the nuclear equation of state for 
supernova simulations and neutron stars. We discuss an updated and improved Notre-Dame- 
Livermore Equation of State (NDL EoS) for use in supernovae simulations. This Eos contains 
many updates. Among them are the effects of 3- body nuclear forces at high densities and the 
possible transition to a QCD chiral and/or super-conducting color phase at densities. We also 
consider the neutron star equation of state and neutrino transport in the presence of strong 
magnetic fields. We study a new quantum hadrodynamic (QHD) equation of state for neutron 
stars (with and without hyperons) in the presence of strong magnetic fields. The parameters are 
constrained by deduced masses and radii. The calculated adiabatic index for these magnetized 
neutron stars exhibit rapid changes with density. This may provide a mechanism for star-quakes 
and flares in magnetars. We also investigate the strong magnetic field effects on the moments 
of inertia and spin down of neutron stars. The change of the moment of inertia associated 
with emitted magnetic flares is shown to match well with observed glitches in some magnetars. 
We also discuss a perturbative calculation of neutrino scattering and absorption in hot and 
dense hyperonic neutron-star matter in the presence of a strong magnetic field. The absorption 
cross-sections show a remarkable angular dependence in that the neutrino absorption strength is 
reduced in a direction parallel to the magnetic field and enhanced in the opposite direction. The 
pulsar kick velocities associated with this asymmetry comparable to observed pulsar velocities 
and may affect the early spin down rate of proto-neutron star magnetars with a toroidal field 
configuration. 



1. Introduction 

To describe the structure and hydrodynamics of compact matter; an equation of state (EoS) is 
needed to relate the physics of the state variables [1]. In supernovae the EoS determines the 
dynamics of the collapse and the outgoing shock, and determines whether the remnant ends up 
as a neutron star or a black hole. In a neutron star, it determines the mass-radius relationship, 
stellar composition, cool-down time and dynamics of neutron star spin down and mergers. In this 



paper we summarize some progress in the development of Equations of state for supernova and 
neutron star simulations. In particular we highlight the role of 3-body forces, QCD, magnetic 
fields and neutrino transport. 



2. Three body forces and the EoS 

At present, only a few hadronic EoSs are commonly employed that cover large enough ranges in 
density, temperature and electron fraction to be of use in core-collapse supernova simulations. 
The two most employed in astrophysical simulations are the EoS of Lattimer & Swesty (LS91) [2] 
and that of H. Shen et. al. (Shen98) OH]. The former utilizes a non relativistic parameterization 
of nuclear interactions in which nuclei are treated as a compressible liquid drop including 
surface effects. The latter is based upon a Relativistic Mean Field (RMF) model using the 
TMl parameter set in which nuclei are calculated in a Thomas-Fermi approximation. Baryonic 
matter was parameterized with a new RMF model that treated nuclei and non-uniform matter 
with the statistical model of Hempel et. al. [5]. 

Here, we discuss a new Notre Dame-Livermore (NDL) EoS This EoS evolves from 
the original Livermore formulation [71 |8], but unlike the previous version this NDL EoS, is 
consistent with known experimental nuclear matter constraints and recent [Hj mass and radii 
measurements of neutron stars. Below nuclear matter density, the conditions for nuclear 
statistical equilibrium (NSE) are achieved at a temperature of T ~ 0.5 MeV. Below this 
temperature the nuclear matter is a approximated by a nine element reaction network which 
must be evolved dynamically. Above this temperature, the nuclear constituents are represented 
by free nucleons, alphas and a single "representative" heavy nucleus. Among the new features 
in the NDL EoS is that high density phase of the EoS is treated with a parameterized Skyrme 
energy density functional that utilizes a modified zero range 3-body interaction. The effects of 
pions on the state variables at high densities is also included as well as the consequences of a 
phase transition to a QGP. 

Above nuclear matter saturation density we include both 2-body {vlj ) and 3-body (v^jl.) 
interactions in the many-nucleon system. The Skyrme two-body potential is given in the 
standard form [9]. 
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Here [B] we consider the possibility that the Skyrme potential can be dominated by a 3-body 
repulsive interaction at high density. This term is taken to be a zero range force of the form 
^^123 = tsS (ri — r2)<5 (ra — rs). If the assumption is made that the neutron-star medium is 
spin-saturated [TU] , the three-body term becomes a density dependent two-body interaction [S] 
that we generalize to a modified Skyrme interaction that replaces the linear dependence on the 
density with a power-law index a. [11] 



vi,' = (l + Ps) S (n - ra) p'^ ("^4^) • (3) 
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A value of o" = 1/3 is a common choice |12J|13J. However, in the present approach we treat o" as a 
free parameter constrained by the skewness coefficient and observed neutron-star properties |14] . 

All quantities and coefficients for symmetric nuclear matter are obtained from the usual 
relations. The pressure, is P = (^)' '^^^ volume compressibility of symmetric nuclear 

matter is calculated from the derivative of the pressure: K = 9 (§^) = 18^ -|- 9n^J^ (^)' 
The skewness coefficient, is from the third derivative of the free energy per nucleon Qq = 



27n^-^ (^) • ^PPlyi'^g the saturation condition P(n = no) = n^-^ \n=no = 0, one obtains 
[6] a system of four equations in terms of Iq, (3ti + 5t2), t^, and a. Solving this system for a 
yields 
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where the subscript zero denotes values at the saturation density. 

For our purposes, we adopt inferred values of uq, Eq, Kq,Qq from the literature and use 
these to determine the Skyrme model parameters. We also demand that these parameters 
allow neutron star masses > 1.97 it 0.04 Mq p!3]. The saturation density uq ~ 0.16 fm~^ 
and the binding energy per nucleon £"0 = —16 MeV are reasonably well established ^15j. The 
determination of the compressibility parameter from experimental data on the giant monopole 
resonance on finite nuclei has been a long standing conundrum. For our purposes we adopt the 
median value and uncertainty from Ref. [16] , i.e. Kq = 240 it 10 MeV as this is appropriate for 
the Skyrme force approach employed here. Solving the saturation conditions self consistently, 
we therefore determine the best range for the nuclear compressibility consistent with the results 

of dE]. 

There is even more uncertainty in the skewness parameter Qo- Breathing mode data [17] 
implies Qo = —700 it 500 MeV. Using the range for Kq given by |16j and solving the saturation 
conditions, we find [6] a skewness coefficient of Qo = —390 it 90 MeV consistent within the 
range given in Ref. [16]. The fiducial NDL EoS is then constructed [6] using the The Skyrme 
coefficients to = -1718 MeVfm^, (8*1 + 5*2) = -102 MeVfm^ = 13226 MeVfm3'^+^ and 
a = 0.369. 

The density dependence of the symmetry energy beyond saturation is highly uncertain. For 
many Skyrme models the symmetry energy either saturates at high densities, or in the worst 
case becomes negative. This results in a negative pressure deep inside the neutron star core. We 
implemented [6] a linearly increasing function of density. The symmetry energy at saturation 
was determined by the difference between the energy per particle for pure neutron matter and 
that of symmetric matter at T = MeV. For all relevant parameter sets the NDL EoS symmetry 
energy at saturation is 5*0 = 30.5 MeV [6]. 



3. QCD and the EoS 

For sufficiently high densities and/or temperature a transition from hadronic matter to quark- 
gluon plasma (QGP) can occur |18) . Progress [19] in lattice gauge theory (LGT) has shown that 
at high temperature and low density a deconfinement and chiral symmetry restoration occur 
simultaneously. In particular, it has been found |19| that the order parameters for deconfinement 
and chiral symmetry restoration changes abruptly for temperatures of T = 145—170 MeV |20tl21j 
as a smooth crossover. 

At low density the hadron phase can be approximated as a pion- nucleon gas, while the 
QGP phase can be approximated in a bag model as a non-interacting relativistic gas of 
quarks and gluons |22j . The LGT results then imply a range for the QCD vacuum energy 
of 165 < S^/^ < 225 MeV. Also requiring that the maximum mass of a neutron star exceed 
1.97 ± 0.04 Mq |14j implies a value for B^^'^ near the top end of that range. 

For the description of quark matter we utilize a bag model with 2-loop corrections, and 
construct the EoS from a phase-space integral representation over scattering amplitudes. We 
allow for the possibility of a coexistence mixed phase in a 1st order transition or a simple cross 
over transition. It is convenient to compute the QGP in terms of the grand potential, 0,{T, V, /u), 
where the grand potential for the quark-gluon plasma takes the form: 



n = + ^92) + ^90 + ^g2 + BV. 
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Where qq and go denote the 0* -order bag model thermodynamic potentials for quarks and 
gluons, respectively, and q2 and g2 denote the 2-loop corrections. In most calculations sufficient 
accuracy is obtained by using fixed current algebra masses (e.g. m„ ~ ~ GeV, 
rus ~ 0.1 — 0.3 GeV). For this work we chose the strange quark mass to be nis = 150 MeV 
and a bag constant B^^'^ = 165 — 220 MeV. The quark contribution to the thermodynamic 
potential is given [18] in terms of a sum of the ideal gas contribution plus a two loop correction 
from phase-space integrals over Feynman amplitudes [23 1. 

Fig. [l] compares [6j the neutron star mass radius relation for the NDL EoS for: 1) a hadronic 
EoS with 3-body forces (solid line); 2) a first order QCD transition with B^^^ = 220 MeV; 
and 3) a simple QCD cross over transition. Also, shown for comparison are results from the 
LS180 EoS, Shen EoS and the original Bowers & Wilson EoS. Note, that all three versions of 
the NDL EoS easily accommodate a maximum neutron star mass > 1.97 it 0.04 Mq, however, 
the hadronic version must have the 3-body forces at high baryon density. A first order phase 
transition to a QGP is consistent with the high maximum neutron star mass constraint jl4] for a 
bag constant B^^^ > 220 MeV. This imposes a low baryon density transition temperature of Tc 
= 158 MeV |22j which is consistent with the current range of crossover temperatures determined 
from LGT \W\. 




Figure 1. Mass-radius relations for the EoS of Shen [3j, Lattimer &: Swesty (LS180) [2], 
Bowers & Wilson [7], compared with the new NDL EoS [6] with and without a mixed phase 
transition or a simple crossover transition to a QGP. Note that the soft LS180 hadronic EoS 
and the previous Bowers &: Wilson EoS without 3-body forces cannot satisfy the astrophysical 
constraint of a maximum neutron star mass > 1.97 it 0.04 Mq |14j as shown. 



4. Magnetic Fields and The Eos 

We have also considered the neutron star equation of state and neutrino transport in the 
presence of strong magnetic fields [23]-f3U]. Indeed, magnetic fields are everywhere in Nature 
and frequently play a role in astrophysical phenomena. In particular the existence of magnetars 
and magnetar flares |3H [32l |33] , along with the observed asymmetry in supernova explosions 
and the observed pulsar kick velocities all suggest the strong magnetic fields play an important 
role in supernova explosions and the formation of proto-neutron stars \34:\ [35l I36j . In view of 
this we have undertaken studies of a variety of phenomena. 

In |25j we considered the role that strong interior magnetic fields {B 10^^ G) would have 
on neutron star structure and stability. We considered the nuclear equation of state for an 
ideal npe gas in a strong magnetic field. In particular, we calculated the proton concentration, 
the threshold densities for neutron, muon, and pion production and pion condensation in a 
strong magnetic field both without and with the effect of the nucleon anomalous magnetic 
moments. It was shown [25j that the higher Landau levels are significant at high density 
in spite of the existence of a very strong magnetic field. In particular, at high density, the 
proton concentration approaches the nonmagnetic limit. In particular, we have obtained the 
neutron appearance threshold density in a magnetic field when the nucleon anomalous magnetic 
moment is included. We also have shown [25j that the muon and pion threshold densities are 
not affected by magnetic fields for B < 10^^ G. We also obtained an equation of state for a pion 
condensate in strong magnetic fields. We found [25] that magnetic fields reduce the amount of 
pion condensation. However, we still could find distinguishable effects from a pion condensate in 
strongly magnetized neutron stars. In addition, we demonstrated an oscillatory behavior of the 
adiabatic index in both strongly magnetized n,p,e and n,p.e,;U,7r gases at high density. Here 
we speculated that this behavior might lead to an interior pulsational instability. 

In [26j we investigated the possibility that soft gamma-ray repeaters (SGRs) and anomalous 
X-ray pulsars (AXPs) might be observational evidence for a magnetic phase separation in 
magnetars. We studied such magnetic domain formation as a new mechanism for SGRs and 
AXPs in which magnetar matter separates into phases containing different flux densities. We 
identified the parameter space in matter density and magnetic field strength at which there is 
an instability for magnetic domain formation. We showed that such instabilities will likely occur 
in the deep outer crust for the magnetic Baym, Pethick, and Sutherland (BPS) model and in 
the inner crust and core for magnetars described in relativistic Hartree theory. Moreover, we 
estimated that the energy released by the onset of this instability is comparable with the energy 
emitted by SGRs. 

In [275 this has recently been extend to a study a new quantum hadrodynamic (QHD) equation 
of state for neutron stars (with and without hyperons) in the presence of strong magnetic fields. 
The parameters were constrained by the condition that deduced neutron star masses and radii 
that must be consistent with the recent observations [U] of a high mass neutron star. The 
calculated adiabatic index for these magnetized neutron stars exhibited the same rapid changes 
with density. This was hypothesized to provide possible insight into the mechanism of star- 
quakes and flares in magnetars. We also investigated the strong magnetic field effects on the 
moments of inertia of neutron stars. The change of the moments of inertia associated with 
emitted magnetic flares was shown to match well with observed glitches in some magnetars. 

In [291 EH] we explored a perturbative calculation of neutrino scattering and absorption in hot 
and dense hyperonic neutron-star matter in the presence of a strong magnetic field. We found 
that the absorption cross-sections show a remarkable angular dependence in that the neutrino 
absorption strength is reduced in a direction parallel to the magnetic field and enhanced in 
the opposite direction. This asymmetry in the neutrino absorption can be as much as 2 % of 
the entire neutrino momentum for an large interior magnetic field. We estimate the associated 
pulsar kick velocities associated with this asymmetry in a fully relativistic mean-field theory 



formulation and show that the kick velocities are comparable to observed pulsar velocities. In 
[3U| we have extended this calculation to include a toroidal magnetic field configuration. In this 
case, there can be an asymmetric emission of neutrino momentum along the magnetic field lines 
that are in the direction of the neutron star spin. This can substantially accelerate the spin 
down of a neutron star in the early cooling phase, ~ 10 sec after core bounce. This is to be 
compared with |28| in which we considered the spin down of a neutron star purely from the 
outflow of neutrinos without a magnetic fleld. 
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